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ABSTRACT 


Equations  are  given  for  ilie  doppler  sliifl  and  rale  of  cliange 
of  doppler  sliift  for  tlie  bistalic  case  wliere  an  orbiting,  nonlrans- 
inllting  earth  satellite  is  illuminated  by  a  transmitter,  and  the 
reflected  energj*  is  received  at  different  locations  on  the  surface 
of  the  earth.  These  equations  have  been  programmed  for  compu¬ 
lation  by  the  N.AREC  computer  for  any  satellite  for  which  the 
orbital  elements  are  known.  The  results  for  a  number  of  satel¬ 
lites  have  been  computed,  using  ir.ansmiitlng  and  receiving  sites 
of  the  Space  Surveillance  System.  Plots  of  various  relationships 
between  doppler  shift,  rate  of  change  of  doppler  shift,  satellite 
height,  earth-center  angle  between  the  receiver  and  the  satellite, 
and  zenith  angle  from  receiver  to  .satellite  are  shown  for  a  typical 
satellite,  1938  Alpha,  E.xplorer  I. 


PROBLEM  STATUS 

This  is  an  interim  reijort  on  one  phase  of  the  problem:  work 
is  continuing. 
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EQUATIONS  FOR  DISTATIC  DOPPLER 
SHIFT  AND  RATE  OF  CHANGE  OF  DOPPLER  SHIFT 
OF  DARK  S.ATELL1TE  OBSERVATIONS 


INTRODUCTION 

In  order  to  separate  reflected  satellite  sipials  from  those  caused  by  aircraft,  meteor 
trails,  Hjtiuninj:  discliarites,  and  direct  transmitter  fcedtiirougii,  the  addition  of  doppler 
sliift  and  rate  of  cltanpe  of  doppler  shift  data  could  he  of  considerable  aid.  Equations  for 
tiie  dopi>ler  shift  f.j  and  rate  of  cltati^e  of  dopitler  sliift  will  now  be  develoitcd  for  the 
histatic  case  where  an  orhitin^  satellite  is  illuminated  by  a  transmitter  at  one  location 
on  the  surface  of  the  earth  and  the  reflected  energy  Is  received  at  a  second  location  on 
the  surface  of  the  earth.  In  particular,  the  transmitter  and  receiver  sites  lie  along  a 
great  circle  which  also  contains  several  other  transmitter  and  receiver  sites.  Tills  com- 
iilc.s  of  transmitters  and  receivers  is  known  as  the  Spaee  Surveillance  System.*' 

The  antenna  configurations  at  each  site  result  in  a  fan-shaped  heam  pattern  with  Its 
wide  dimension  in  t.he  jilane  of  the  great  circle.  The  heam  jiatterns  overlap  such  that 
whenever  a  satellite  crosses  the  great  circle  jilane  it  must  iiass  througli  one  or  more 
heam  patterns.  This  results  In  a  -fence'  of  beam  patterns,  and  for  this  reason  the  great 
circle  |)lane  is  referred  to  as  the  fence  jilane.  .At  each  receiver,  data  is  recorded  and 
analyzed  to  determine  time  of  jiassage  of  the  satellite  through  the  fence,  the  zenith  angle 
of  the  satellite  from  the  receiver  at  the  time  of  the  observation,  and  identity  of  the 
satellite. 

It  more  than  one  receiving  st.ation  simult.ineously  observes  the  same  satellite  jias- 
sage,  its  height  may  be  determined  by  triangulation. 

The  coordinate  system  used  in  tlie  derivation  of  the  equations  is  defined  in  Fig.  1. 

The  earth  is  oriented  such  that  the  X  a.sis  Is  along  the  longitude  of  Greenwich  (O'*),  the 
Y  a.\is  is  along  90''W  longitude,  and  the  7.  a.sis  is  along  the  north  polar  a.\is.  Thus,  the 
\V  Jilane  corresjionds  to  the’earth's  equatorial  jilane.  Longitudes  are  measured  jiosilive 
west  ot  Greenwich,  .md  latitudes  are  measured  jiositive  towards  the  7  a.sls.  Tlie  coor¬ 
dinate  a.ses  are  fi.scd  with  resjiect  to  the  earth  and  must  therefore  rotate  with  ilic  earth. 

Figure  2  shows  the  orientation  of  the  fence  plane  and  its  great  circle  Intersection 
with  the  surface  of  the  earth.  The  locations  of  a  receiver  R  and  transmitter  T  along  the 
great  circle  are  also  shown. 

The  geometry  of  a  satellite  observation  in  the  jilane  of  the  fence  is  shown  in  Fig.  3. 
The  sejiaration  of  transmitter  and  receiver  is  e.xaggerated  for  clarity. 


PROCEDURE 

The  .NAREC  computer  was  programmed  to  compute  the  doppler  shift  f,,  and  the  rale 
of  change  of  dojijiler  shift  f,,  for  any  satellite  for  which  the  orbital  elements  are  given. 
The  fence  plane  may  be  any  jilane  for  which  the  equation  is  given,  and  the  receiver  and 


‘i'Proc.  of  the  IRE  ■ISf.N'o,  •}):663-669,  “The  Navy  Space  Surveillance  System,"  R.  L.  Easton 
and  J .  J .  Fleming. 
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irniisiniiicr  may  be  locniccl  anyvvlicrc  along  ilic  fence  great  circle.  For  all  the  answers 
com|)utc(l  to  date,  the  fence  plane  Is  that  of  the  Space  Sureeillnncc  System.  The  receiver 
is  located  at  either  Ft.  Stewart,  Georgia,  or  Silver  Lake,  Mississippi,  ami  the  transmitter 
is  located  at  Jordon  Lake,  Alabama.  Tiie  sulxiivision  of  the  region  above  the  transmitter 
and  receiver  is  illustrated  in  Fig.  4,  using  the  data  from  the  satellite  1958  Alpha  as  an 
e.Nample.  Again,  the  separation  between  the  receiver  and  transmitter  is  e.Naggerated  for 
clarity. 


Fig,  •!  -  Sub(!:v:s:o!'.  o:  she  region 
above  the  sr.'i.nsrsiitter  ar.d  receiver 


The  first  compulation  was  for  the  case  whore  the  satellite  Is  located  at  iJolnt  1  at  llic 
lime  It  crosses  the  fence  plane.  At  point  1,  the  starting  height  ll„  Is  220  statute  miles 
and  the  starting  earih-conicr  angle  is  T^W.  At  the  com|)lciion^of  the  com|)uiatlon  of 
fci,  and  f,i|  at  |)oiin  1,  att  increment  of  .ill  =  50  miles  is  added  to  ll^  and  fn  and  f.i,  for 
point  2  are  com|)uiod.  Another  increment  of  .:ii  Is  then  added  and  the  process  repeated 
until  a  finishing  height  of  11,^=  1220  miles  is  reached  and  f.i,,  and  are  computed. 

At  this  time,  the  height  is  returned  to  ,  an  increment  of  =  l”  is  added  to  ,  and  f,, 
and  f,!,..  are  computed.  In  this  manner,  ! ,  and  f,,  are  computed  for  each  of  Uio  points  '  ” 
1  through  231. 

hi  the  usual  application  of  the  program  to  the  computation  of  f  and  f .  for  actual 
saicliitos,  the  value  oif  ll„^  is  chosen  as  the  height  of  perigee  and  ll’„  is  chosen  as  the 
height  of  apogee.  Fo^r  all  computations  made  to  date,  has  boon  chosen  as  7%  of  the 

iocei\or  and  Vo^  =  3  L,  although  any  value  of  and  may  bo  used.  Incronionis  of  .111 

and  are  chosen  to  obtain  the  desired  subdivision  of  the  region  of  interest,  hi  the  case 
iliusiratod,  there  are  21  incremonis  of  height  and  11  iiicrcmcnts  of  oarih-conlor  angle 
giving  a  231 -point  coverage  of  the  region.  ’ 

For  each  point  al  which  a  compulation  is  performed,  four  cases  must  bo  considered 
depending  upon  whollier  the  fence  crossing  is  in  a  north -to-soulh  or  south -lo-norlh  dircc- 
uon,  and  upon  wholhor  tlio  saiollile  is  approaching  apogee  or  perigee  in  its  orbit: 

Case  1.  North -south,  approaching  apogee  (NSAA); 

Case  2.  North -south,  approaching  perigee  (NSAP); 

Case  3.  South-north,  approaching  apogee  (SNAA); 

Case  4.  South-north,  approaching  perigee  (SNAP). 


The  program  is  written  to  compute  answers  for  Case  1  first,  then  to  return  to  point  1 
and  compute  answers  for  Case  2.  hi  a  similar  manner,  answers  for  Cases  3  and  4  are 
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then  computed.  If  predictions  are  given  for  the  satellite  of  interest,  the  appropriate  case 
will  be  known.  Otherwise,  there  are  four  sets  of  f.,  and  (,  answers  for  each  point. 

After  the  computer  finishes  computing  the  answers  for  Case  4,  the  computer  stops 
and  a  new  data  tape  is  read  into  the  computer  with  the  elements  of  the  ne.xt  satellite  for 
which  answers  are  desired. 


CO.MPUTER  l.N'PUTS  A.\D  OUTPUTS 


Data  tapes  have  been  insertetl  and  answers  obtainetl  for  the  following  lists  of  satellites: 
A.  With  the  receiver  located  at  Ft.  Stewart  and  the  tr.msmltter  at  Jordon  Lake: 


1.  1958  Alpha 

2.  1938  Beta  2 

3.  1939  Alpha  I 
•1.  1939  Epsilon  2 
3.  1939  Eta 

6.  1939  Iota  1 

7.  19G0  Beta  2 

8.  Circular  Orbit  with  II  = 


E.xplorer  I 
Vanguard  I 
Vanguard  II 
Discoverer  V  Capsule 
Vanguard  III 
E.Nplorer  Vll 
Tiros  1 

100  miles:  inclination  =  90®. 


B.  WiUi  the  receiver  iocated  at  Silver  Lake  and  the  transmitter  at  Jordon  Lake: 


1. 

1938  Alpha 

E.xplorer  1 

2. 

1939  Alpha  1 

Vanguard  II 

3. 

1939  Epsiion  2 

Discoverer  V  Capsule 

4. 

1939  lota  1 

E.xplorer  Vll 

3. 

1960  Beta  2 

Tiros  I 

6. 

1960  Gamma  1 

Transit  IB,  Second  Stage 

7. 

1960  Epsiion  1 

Sputnik  IV 

8. 

1960  Zeta  1 

Midas  11 

9. 

1960  Eta  3 

Transit  IIA,  Second  Stage 

10. 

1960  Iota  1 

Echo  I 

Figure  5  is  the 
_ _  _ 1  .  . 

first  page  of  the  printout  of  the  results  for  1958  Alpha 

- - — —.  The  answers  are  divided  into  several 

groups,  two  of  which  arc  shown.  The  first  group  includes  all  the  answers  for  each  of  the 
21  height  increments  for  -/=  7  W.  The  second  group  includes  the  answers  for  the  21 
height  increments  for  y=  6  W.  Tlie  succeeding  pages  of  answers  (not  shown)  give  answers 
for  each  of  the  other  increments  of  earth-center  ajigie.  The  positions  of  the  decimal 
points  in  the  answers  are  shown  in  the  first  group,  and  are  in  the  same  positions  in  all 
groups.  The  identity  of  the  answers  in  the  printout  may  bo  determined  by  referring  to 
Fig.  6. 
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uherc  i  «  timv,  in  seconds,  Tor  which  all  equations  are  evaluated  (t  »  Ofor  all 
contiiutations  made  to  date): 

I  =  iKcltnation,  In  dettrees,  of  the  orbital  plane; 

'  =  eccentricity  of  the  orbit: 

»  =  semtmajor  axis,  in  miles. of  the  orbit: 

•  s  earth-center  anyle  lietween  receiver  and  satellite,  “W  -♦ 

L„=  latitude  of  satellite  at  first  height  In  the  group: 

„  =  longitude  of  satellite  at  first  height  In  the  group; 

=  zenith  xtglc,  in  degrees,  from  receiver  to  satellite; 

'•'sF  ■  distance,  in  miles,  between  receiver  and  satellite; 

'*si  ■  distance,  in  miles,  between  transmitter  and  satellite; 

"  rate  of  change  of  radius  vector.  In  miles  per  second  (also  equal  to 
rate  of  change  of  height): 

Lj  =  rate  of  change  of  latitude.  In  degrees  per  second: 

“s  ■  change  of  longitude.  In  degrees  per  second; 

f,i  =  dopplcr  shift,  in  cycles  per  second: 

f,,  =  rate  of  change  of  dopplcr  shift,  in  cps  per  second. 

Fig.  6  -  Code  of  the  answers  on  the  computer  printout  (Fig.  5) 
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Certain  i:iformatioii,  :sucJi  as  orljilal  cleincmf!,  starting;  and  flinsluiiu  heights,  and 
starting  and  ftnisliing  earth-center  angles,  is  required  lor  ilie  data  tape  iniiui  to  the  cont- 
puier  for  eacli  satellite  for  which  answers  are  desired;  the  inforinaiion  required  is  ll:o 
following: 

e,  eccentricity  oi  the  orbit: 

a,  seiniinajor  axis  of  the  orbit,  in  statute  miles: 

i ,  inclination  of  the  orl)ital  iiiane,  in  degrees: 

1,  time  for  which  answers  are  to  he  cuni|)uled; 

»,  seiccted  toierance  within  which  successive  iterations  in  the  appruximalion  to 
the  soiution  to  Kepler's  equation  must  fall: 

■'os>  starting  earth-center  angle,  in  degrees: 

>0,1  finishing  earth-center  angle,  in  degrees; 

,  seiccted  increments  of  earth-center  angle.  In  degrees: 

l!„^,  starting  height,  in  statute  miles; 

ll^,,  finishing  height,  in  statute  mites: 

.'.11,  selected  Increments  of  height,  in  statute  miles. 

hi  addition,  certain  other  data  are  included  on  a  second  input  tape  which  contains 
quantities  which  are  either  alw.ays  consiam  or  remain  constant  during  several  changes 
of  the  data  tapes.  The  information  on  the  Constants  Tape  Includes  the  following: 

1.^,  latitude  of  the  receiver,  in  degrees; 
longitude  of  tlic  receiver,  in  degrees: 

l.j,  latitude  of  the  transmitter,  in  degrees: 

T,  longitude  of  the  transmitter,  in  degrees: 

A,  coefficient  of  the  \  term  in  the  equation  of  the  fence  plane: 

It,  coefficient  of  the  y  term  in  the  equation  of  the  fence  plane: 

C,  coefficient  of  the  t  term  in  the  equation  of  the  fence  plane; 

radius  of  the  earth,  in  .siainic  miles; 

angular  rotation  rate  of  the  earth,  in  radians  iier  mean  solar  second: 

”,  the  constant. 

RESULTS 

Some  of  the  results  for  1958  Alpha,  Case  1  (norlli-to-south,  approaching  apogeob 
ai  e  shown  graphically  in  Figs.  7,  8,  and  9.  Figure  7  is  a  plot  of  rate  of  change  of  doppler 

shift  vs  doppler  shift  for  a  family  of  curves  of  constant  height  and  a  family  of  curves  of 

constant  eaith-contor  angle.  The  data  for  Case  2  (north-to-soutli,  approaching  perigee), 
if  plotted  on  the  same  graph,  would  result  in  a  closed  curve  for  each  curve  in  the  family 
of  constant  earth-center  angle.  For  clarity,  only  the  data  for  Case  1  are  shown. 

Figure  8  is  a  plot  of  rate  of  change  of  doppler  shift  vs  zenith  angle  for  a  family  of 
curves  of  constant  earth-center  angle.  As  with  Fig.  7,  the  data  for  Case  2  could  also  be 
included  on  Fig.  8  to  form  a  family  of  closed  curves. 
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E'ig.  ”  -  K.itf  o:‘  change  of  doppicr  shift 
shift  for  curves  of  cor.st.ir.t  heigh:  I'l 
e.irth-cef.ter  atjgJe 


vs  (ioppier 
;«i  con  star.: 


Fig.  8  -  Rale  of  change  of  doppler  shift  vs  zenith  angle  for 
curves  of  constant  height  and  constant  earth-center  angle 
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F:g.  V  -  Height  vs  zer.t'.h  artglt  /or  curves 
o:  const.if.t  tioppier  sh::t 


By  means  of  Figs.  7  and  8,  data  were  obtained  to  plot  Fig.  9,  a  plot  of  lioiglil  vs  zoiillli 
angle  for  a  family  of  curves  of  constant  f,,. 

The  computer  i)rintouts of  answers,  one  page  of  which  is  shown  in  Fig.  5,  have  boon 
comi)iled  into  two  convenient  booklets  {an  original  and  a  carbon  coity)  for  each  of  tlie 
eighteen  comimter  runs.  The  answer  tapes  are  on  file  from  which  any  number  of  print¬ 
outs  may  be  made  using  any  Fle.xowriter  with  a  carriage  105  or  more  S|)aces  wide. 


EQU.^TIO.NS 

The  equations  wliicli  were  used  in  the  comi)uter  program  will  now  be  given.  Symbols 
used  in  the  equations  and  in  the  te.xt  are  summarized  in  Ai)pendi.x  A.  The  doppler  shift 
and  rate  of  cliange  of  doi)i)ler  shift  are: 


(Jopp!  cr  shift: 


'ST' 


rntc  of  rli.Tngo  of  doppler  shift: 


+  H, 


ST 
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wliore 


=  iranstnitlcd  fre(|uoiicy: 

<'=  prupatLation  cunsiaiil: 

•ij.,.  =  (iistaiice  Ijeltt-ecMi  saiclliic  and  receiver: 

<!^^=  (iistaiice  lieltt-een  satellite  anti  transmitter. 

Time  varialile  expressions  for  <'};,{  and  <*5;^  ‘'•’e  recpiired  so  that  they  may  he  differ 
enliati'd  to  ohtaiii  their  first  .and  second  time  derivatives: 


and 


•  V'  ■  -V*’  ■  >B-‘  ■  '  "r  ■)'  ' 

j-V*)  ■  ■  iV)  ■  'i!'  ■  V')  ■ 


’.vliere  j  Vjj(:),  Vj.(!),  are  the  time  variable  rectanitle  coordinates  of  satellite  position: 
l''K'  '■»(’  ^rI  rectan^ilar  coordinates  of  the  receiver:  and  ('t^,  y^,  7^)  are  the 

rectan:4t!lar  coordinates  of  the  tiMiisniitter.  The  rectangular  coordinates  of  satellite  [losl- 
llon  are  obtained  front  the  spherical  coordinatos  by  the  following  transfonnatlons: 


y-)  =  fgd)  -j-n) 


•V'f  '  'sf  '  ■s'''' 


V'>  '  'sf*' 

where  rj.(i)  is  the  time-dependant  radius  vector  of  the  satellite,  and  ( '5(1),  ..5(1),  i-sCt)) 
arc  the  time-dependant  direction  cosines  of  satellite  iiosillon.  Tlie  direction  cosines  of 
satellite  position  are  defined  by  the  following  Uircc  simultaneous  equations: 

n  •  •  ^sfM  I 


•vf)  -sO  •  -s<'>  •  =s('>  =  0 


•;ft)  -.sC) 


..;(!)  -sd)  -  =  cos  .Jft) 


where  ( <j(i)j  are  the  coefficients  of  the  satellite  plane:  [  ■/(i),  -;(«),  ■-'o(«)l 
arc  the  direction  cosines  of  the  satellite  after  the  earth’s  rotation  but  before  o”rbital "motion 
and  ;(i)  is  the  earth-center  angle  through  which  the  satellite  moves  in  its  orbital  plane. 

The  coefficients  in  the  equation  of  the  satellite  plane  are  found  as  follows: 


(A2 

>  152  i 

B'd) 

(A2 

t  ^ 

C'(t) 

(A^  H 

^  n2  + 

NAVAL  RESEARCH  LABORATORY 


11 


Whore  (A'{t),  a‘{t),  c‘{i)l  are  the  coefficieiils  in  the  equation  of  the  satellite  plane  before 
normalization.  Those  coefficient.*^  are: 


A*(l)  -  -It  «in  4_I  •  A  Cos  I 

p  p 

lt*(t)  -  A  sin  t  •  It  cos  i  I 
P  p 


C'M)  ^  C 

Whore  (A,  B,  c)  are  the  coefficients  in  the  equation  of  the  satellite  plane  at  t=  0: 


=  angular  rotation  rate  of  earth, 
i-,-  =  precession  of  ihe  plane  of  the  satellite, 
t  =  elapsed  time  since  fence  crossing, 

A  = 

H  =  -K:  „  ,  :UUi 


where  (  :>„)  are  the  direction  cosines  of  the  sateilite  at  «  =  0,  and  K  is  defined  by 

the  following  quadratic  equation: 


5in* 


i  i)  K*  -  (2-.  ..  sin  i)  K  - 


Sin*  I 


i)  =  0, 


Hero,  i  is  the  inclination  of  the  sateiiite  plane.  Tlie  oartit-center  angie  =  through  which 
the  sateiiite  moves  in  its  orbital  plane  during  time  i  is  given  by  the  equation 


aft)  =  Vg(l)  - 


where  ^5(1 )  =  nme  variable  o.xpression  for  true  anomaiy:  v„=  true  anomaly  at  i  =  0: 
and  !.  ^=  angular  rotation  rate  of  perigee.  Now, 


V'sft) 


nrccos 


rcos  Esfl)  -  c) 
1  -  f  cos  Egf  I  )J 


where  £5(1)  -  time  variabie  e.xpression  for  eccentric  anomaly,  and  e  =  eccentricity  of  Uie 
orbit.  To  obtain  £5(1),  the  transcendental  equation  known  as  Kepler’s  equation  must  be 
solved: 


Es(‘)  =  Msf')  *  S'" 


where 


Ms(t)  =  n{t  -  t)  =  time  variable  expression  for  mean  anomaly, 

2n 

n  =  y  =  mean  angular  velocity  of  sateiiite  in  its  orbit. 


P  =  anomalistic  period  of  the  orbit. 
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''  ”  ■  ~  =  f*'*'*-*  ‘>1^  perigee, 

“  *o  ■  '■  =  »>eaii  aiiumnly  at  i  =  0, 


i:  = 


eccentric  anomaly  at  «  =  0, 


a  =  seininiajor  axis  of  the  orbit,  and 
=  radins  vector  of  tlie  satellite  at  i  =  0. 

Hio  direction  cosines  of  tlio  satellite,  considering  the  earth's  rotation,  but  before  orbital 
motion,  are  defined  as  follows: 

•.;f  «)  =  ro^  co«  t) 

t)  =  ens  l.^r  t)  <in 
:  ;f«>  =  tin  l.;f  t) 


wlieri 


Also, 


=  angular  rotation  rate  of  earth,  and 
.  =  precession  of  the  plajie  of  the  satellite. 

I,  =  nrcsin  :  ,  nnd 


*•  -  nrcsin 


i„2  -  ■•2)  '■'V 

r  o  o'  / 


where  (  are  Uie  direction  cosines  of  the  satellite  at  t  =  0.  The  direction  cosines 

‘  'o’  ••o>  'o  '  satellite  at  t  =  0  ore  defined  by  the  following  three  simultaneous 

eciiiaiions: 

>.2  •  u2  -  1.2  =  I 
o  o  o 

‘’f  '-o  *  '’f  ''o  ^  ''f  ‘'o  =  0 

^  ^  "o  "  ■'r  ’-"o  =  '/o 

where  [hj.,  b^,  c^.]  are  the  coefficients  of  the  fence  plane  equation;  /i  ^>  ]  are  the 
direction  cosines  of  the  receiver;  and  is  the  earth-center  angle  between  receiver  and 
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satellite  at  «  =  0.  Given  the  satellite  height  ii„  anti  ihe  senitli  angle  7.^  from  the  receiver, 
the  earth-center  angle  may  be  found  by  means  oi  the  equation 

Iv  •‘h* 

1.*,-  II 

where  is  the  earth's  radius.  Hie  time  variable  e.xpression  for  satellite  hciglit  is: 

iigdi  -  rs(«)  -  R, 

where  r^ti)  &  a|  t  -  v  cfl<i  is  the  radius  vector  of  the  satellite,  »and  v  arc  as  defined 
before,  and  ~  time  variable  e.\|)ression  for  eccentric  atiomaly  (previously  defined  by 
Kepler's  equation). 


CONCLUSIONS 

As  Implied  by  the  plot  shown  in  Fig.  7,  doppler  sliift  ajtd  rate  of  change  of  doppler 
shift  are  sufficient  to  define  the  imsltion  of  a  satellite  for  which  semlmajor  a.\ls.  Inclina¬ 
tion,  and  eccentricity  are  known.  However,  due  to  the  relatively  short  duration  of  time 
that  a  satellite  Is  normally  In  Uie  beam  of  the  Space  Surveillance  System,  It  Is  difficult 
to  obtain  a  sufficiently  accurate  measure  of  the  rate  of  ehange  of  do|)pler  shift  to  make  use 
of  a  plot  of  tills  kind.  For  other  systems,  where  several  seconds  or  more  arc  available, 
rate  of  change  of  doppler  shift  could  be  measured  with  sufficient  accuracy  to  define 
satellite  imsltlon. 

Since  the  zenith  angles  of  obseri-ntlons  are  determined  from  other  data  recorded  by 
the  Space  Surveillance  System,  the  addition  of  dop|)ler  information  to  this  data  makes 
|)osslble  Uie  determination  of  the  heigiit  of  a  satellite  with  known  orbital  elements  from 
an  observation  by  a  single  station,  rather  than  requiring  trlajigulatlon  by  coincident 
observations  by  two  or  more  stations.  A  plot  such  as  that  shown  In  Fig.  9  would  be  use¬ 
ful  for  this  purpose.  .A  program  modification  could  be  made  such  that  the  height  corre- 
S|X>nding  to  an  observed  zenith  angle  and  doppler  shift  could  he  computed  directly,  so 
that  It  would  not  bo  necessary  to  refer  to  a  plot.  Otherwise,  four  plots  (one  for  each  of 
the  four  cases)  must  be  prepared  for  each  known  satellite.  This  task  could  be  simplified 
by  having  Uio  NAREC  computer  answer-tape  punched  in  the  proper  format  so  Uiat  the 
plots  could  be  made  on  tiie  automatic  plotter. 

E.xperimental  doppler  data  suitable  to  verify  the  results  have  not  been  available. 
Various  program  modifications,  such  as  indicated  below,  may  be  incorporated  later 
should  such  data  indicate  the  necessity. 

1.  Consider  height  above  mean  earth’s  radius  of  the  receiver  and  transmitter. 

2.  Correct  earth’s  radius  for  oblateness  of  the  earth  when  computing  satellite  heights. 

3.  Consider  decay  of  the  semimajor  axis. 

If  the  coordinates  of  both  tlie  transmitter  and  the  receiver  are  equal,  implying  that 
both  are  at  tlie  same  location,  the  bistatic  case  reduces  to  the  monostatic  case,  for  which 
the  equations  given  in  tliis  report  hold  equally  well. 
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APPKNDIX  A 

DKFINITIOX  OF  SYMllOl^S 


h(.  !»,.  «•,! 

(A.  It.  cj 
iA'(t).  »'(«). 

f;n)! 


V') 

**ST^ 


O 


EsO) 

f,i<0 

ft 

:.ii 


V «) 


Suininiajor  axis  (jifvcii  as  one  of  llio  orbiial  elements) 

Coefficients  of  the  fence  plane  equation  (given) 

Coefficients  of  the  satellite  plane  equation  at  t  -  0 

Coefficients  of  tlie  satellite  plane  equation  considering  the  eartlFs 
lotation  anti  precession  of  ilte  notle 

Nonnallitcti  coefneients  oi  tlic  satellite  plane  equation  consitloring 
tlic  cartli-s  rotation  and  iirecession  of  the  node  ^ 

l»ro|)agatlon  constant  (given) 

IJistanco  hettvecn  satellite  anti  receiver 

Distance  between  satellite  and  transmitter 

Eccentricity  (given  as  one  of  ilie  orbital  elements) 

Eccentric  nnotnnly  at  i  =  0 

Eccentric  anomaiy  at  any  time  t 

Doppler  shift 

Hate  of  change  of  dop|)lor  shift 
Transmitted  frequency 
Selected  Increments  of  height 
Height  of  the  sateilile  at  t  =  0 

Starting  height  of  the  sulxilvlsion  of  the  region  above  the  trans¬ 
mitter  and  the  receiver 

Finishing  height  of  tlie  suixiivision  of  tlie  region  above  liie  traiis- 
mttter  and  the  receiver  uie  u.us 

Height  of  the  satellite  at  any  time  i 

Inclination  of  the  satellite  plane 

Constant,  defined  by  tlie  quadratic  equation: 


sin2  i  -  cos2  i)  -  (2.\„  sin  i)  K  .  (>.2  sin^  i  -  ,.2  cos^  i  )  = 


Lq  Latitude  of  the  satellite  at  i  =  0 

t )  Latitude  of  the  satellite,  considering  only  the  eartli’s  rotation 
and  precession  of  tlie  node 

Lsft)  Latitude  of  the  satellite  at  any  time  t 

Mq  Mean  anomaly  at  t  =  0 

V  t)  Mean  anomaly  at  any  time  t 
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IG 


I) 

I* 


re(i) 


Vj;(  «  J 

I-VsC). 

•  'b*  'k'  *lf 

't'  '"t”  ^T' 
7 

O 

If  n 


\ 


-V’' 


•  K-  -K’  •  K' 

•  *o*  "o* 

-■fn.  -.-(n: 

‘r 


Mean  aiiiiixilar  velocity  of  the  salelliio  in  its  orbit 
Anuinalistic  iioriotl  (uiven  as  one  of  the  orbital  elements) 

Uaclins  vector  of  the  satellite  at  t  =  0 
Hadiiis  vector  of  the  satellite  at  any  time  i 
Uatlius  of  the  earth  ({tiveii) 

True  anomaly  at  •  =  0 
True  anomaly  at  any  time  i 

Hectaimular  coordinates  of  the  satellite  at  any  time  i 

Hectantitilar  coordinates  of  the  receiver 

Uectanitular  coordinates  of  the  transmitter 

Zenith  anj;le  of  a  satellite  as  observed  from  the  receiver  (given) 

Angle  In  tlie  plane  of  the  satellite  through  which  the  satellite 
moves  during  time  i 

Selected  iticreinents  of  earth-center  angle 

Kartb-center  :mgle  between  the  satellite  and  the  receiver  at  t  =  0 

Starting  eartli-cetiter  angle  of  the  snlxlivlslon  of  the  region  above 
the  trausinltter  and  die  receiver 

Finishing  earth-center  angle  of  the  sulxllvision  of  the  region  above 
the  utinsmitter  and  the  receiver 

Tolerance  within  which  successive  Iterations  bt  the  approshnatlon 
to  the  solution  of  Kepler's  equation  must  fall 

Longitude  of  the  satellite  at  t  =  0 

Longitude  of  the  satellite  at  any  time  t  considering  only  the  earth's 
rotation  and  precession  of  the  node 

Longitude  of  the  satellite  at  any  time  t 

Direction  cosines  of  the  receiver 

Direction  cosines  of  the  satellite  at  t  =  0 

Direction  cosines  of  the  satellite  at  any  time  i  considering  the 
earth's  rotation  and  precession  of  the  node 

Direction  cosines  of  the  satellite  at  any  time  \ 

Time  satellite  was  last  at  perigee  before  t  =  0 

Angular  rotation  rate  of  the  earth  (given) 

Apparent  angular  rotation  rate  of  the  satellite  plane 

Precession  of  the  node  (given) 

Rotation  of  perigee  (given) 
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